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One dimensional metals are described by Luttinger liquid theory. Recent experiments have 
£ — ' addressed the relation between this non-Fermi liquid behavior and the existence of a Fermi surface. 

^\ , We show that Luttinger's theorem, with few modifications, holds for the one-dimensional Tomonaga- 

CTN ' Luttinger model. The implications for the high temperature superconductors are discussed. 
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t, PACS numbers: 71.27.+a,74.70.Kn,74.20.Mn 

Soon after the discovery of the high temperature superconductors (HTSC), Anderson jl| noted that the two di- 
mensional (2D) copper-oxide planes in these materials are responsible for the high critical temperatures. He suggested 
that the essential properties of the HTSC's are contained in the physics of the 2D strongly correlated electron liquid. 
While there is a consensus on this question, there are debates on the nature of the low energy physics of that liquid. 

The question is whether 2D quantum fluids are described by Landau's Fermi liquid theory (FLT) or by a theory 
which resembles the physics of one dimensional (ID) systems with short range interactions |^||,[|. Perturbation 
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theory shows that despite some peculiarities compared to three dimensions |6|,0, the essential physics of the 2D 
» , weakly coupled electron liquid is described by Fermi liquid quasiparticles. However, in ID perturbation theory is 
known to violate some of the exact results and therefore it is feasible that in 2D the same could happen leading to 
unknown behavior. ID metals with short range interactions are described by Haldane's Luttinger liquid theory || 
(LLT). In ID, the exactly solvable Hubbard || and Tomonaga-Luttinger [ p~0|Jl~l] ] models capture the essential physics, 
where instead of a simple pole, the one-particle Green's function has a branch-cut singularity and spin and charge 
propagate with different velocities |l2|,|l3| (the Hubbard model in the metallic phase scales (in renormalization group 
sense) to the Tomonaga-Luttinger model |l4] [l6| ) . These properties of ID systems are also important in understanding 
the properties of the organic and inorganic quasi-lD metals. |T^,|l8|. While we know that in 2D when very strong or 
long range interactions are present the Fermi liquid picture breaks and the system can develop charge or/and spin 
■ density waves, or condense into a Wigner crystal, the question is if this is the case for arbitrary small, short range 
\ interactions as is the case in ID. 

Some time ago, Luttinger | fl9f proved a theorem using perturbation theory showing that the volume enclosed by 
the Fermi surface is an interaction strength invariant (therefore the Fermi sphere is incompressible) and that the 
momentum distribution function has a discontinuity at the Fermi momentum of the noninteracting system [^of . In 
ID, there is no discontinuity in the momentum distribution function at, pp, the Fermi momentum of the noninteracting 
system and the excitations consist only of collective boson modes. The momentum distribution function, in the vicinity 
of the Fermi "surface", behaves as |2~lf| \k — pp\ a with nonuniversal, coupling dependent a, and no quasiparticles are 
present in the liquid. The absence of quasiparticles near the Fermi surface gives the distinct properties of the ID 
metals. Nevertheless, even in ID, the Fermi momentum of the noninteracting system plays an important role when 
the interactions are switched on. The exact relationship between Fermi liquid behavior and the existence of a Fermi 
surface is not yet clear and we feel it is important to discuss this relation in Luttinger's theorem. In the case of ID 
metals there have been controversial statements [|2 -25 about the validity of this theorem and this is a question that 
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O . we would like to address in this short paper. 
O ■ In this letter, we prove Luttinger's theorem in ID and therefore show that the theorem holds in a system in which 
the interactions do not produce quasiparticles: the ID g-ology model pl| . This theorem has been proven order by 
order in perturbation theory using the general properties of a Fermi liquid-like Green's function. However, in ID 
there exists a closed integral equation for the single-particle Green's function obtained after the summation of the 
perturbation series. A particular case permits an explicit solution in k — space for the Green's function which is 
a double-valued function of the frequency. No k — space solution is available in the general case, but few results 
regarding the single-particle Green's function can be proven without having an explicit expression for its solutions. 

Ignoring umklapp and the backward scattering processes, the g-ology model describes a set of electrons in ID with 
the Hamiltonian fl2ll]: 



* e-mail blagoev@ph01.bc.edu 
^e-mail kevin.bedell@bc.edu 



1 



a k2,l3 a k2+P,/3 a ki —p,a 



H = Hq + H mt (1) 
H = ^2 v F (k - p F )a\ a ak,a + ^ v F (-k - p F )b\ a b K<7 (2) 

k,(T k.a 

Hint = — [T^a^^bl^pbk^p^ci^-p^ 

+ 6 I 1 ,a 6 L,/3 & fe2+P,/3&fc 1 ^a)] (3) 

where a^,a (6^,6) are the creation and annihilation operators for particles at the +p F {~Pf) branch respectively. 
(3 ~ 9i\\^a,/3 + 9iJ-$a,— f3 (i = 2,4) and || , and _L correspond to particles with parallel and antiparallel spins 
respectively. Using a Ward identity, one obtains Dyson's equation as a singular integral equation for the single 
particle Green's function for electrons moving to the right : 

G+{p,e) = G Q+ {p,e) 

X t 1 + A/ / dkduG + (p-k,e-u)K(k,u)], (4) 



where 



Bi 
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and 



u — Uik + iS[k] oj + Uik — iS[k] J ' 



(5) 



Go+O,e) = / \ , rjr r (6) 

e - v F (p - p F ) + ib[p - p F \ 

is the single particle non-interacting Green's function, and the constants A4 and Bi depend on the couplings g 2 || , g2± 
(for particles on different branches) and 5411 , g4± (for particles on the same branch) [2l| l (we assume that the couplings 
are momentum independent). The charge and spin velocities, u c and u s respectively, are functions of the coupling 
constants. Here p F = Tm/2 and S[q] = Ssign[q]. 



Luttinger's theorem [^9 20 states that in d dimensional space: i. 

2 f dp = 2 f 9(p - p F )dp =V F = (2Tr) d ^ (7) 

JG(p,0)>0 J v 

where p F is the Fermi momentum of the noninteracting system, N is the mean number of particles in the system, V is 
the volume of the system, V F is the volume of the Fermi sphere and G(p, 0) is the interacting Green's function, ii. the 
momentum distribution function u(p) has a discontinuity at the points {pf}> at which the noninteracting distribution 
function n (p) ~ 9(p — p F ) is discontinuous. The discontinuity of n(p) is proportional to the quasiparticle residue, 
i.e. lim [n(p < p F ) — n(p > p F )] = Z. However, in ID Z is zero and the generalized statement is that the derivative 

of n(p) is singular with a power law singularity instead of a delta function singularity. We will show that Luttinger's 
theorem is satisfied (in ID) by the Green's function satisfying Eq.(|J), i.e. that it changes sign when crossing the Fermi 
momentum of the noninteracting system. In this case the interacting Green's function is a product of (see Eq.([|)) the 
noninteracting Green's function, which changes sign at p F and a term which, if it does not change sign when we cross 
p F , Eq.(Q) and therefore the first part of Luttinger's theorem will be satisfied. From now on we work with particles 
on the right branch, i.e. moving to the right. Let us denote by D(p) the term in the parentheses at e = 0, i.e. 

D (P) = ^J J dkdcjG+{p-k,-ij)K(k,u;) (8) 

We can integrate (we explain later in the paper in more detail, how these integrals are calculated) over u using the 
general properties of fermionic Green's function p^| : i. The singularities are located in the second quadrant in the 
complex-frequency plane for p — p F < and in the fourth quadrant for p — p F > 0; ii. G(p, uo) — . The result for 

p > p F is 
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— V [Bi / G+( P -k, Ul k)dk 

—Bi / Gj r {p—k^u i k)9{pF~P J rk)dk 
Jo 

poo 

+Ai / G + (p-k,-u. l k)9{p-k-p F )dk] (9) 



and for p < pp 



G+ (p — k, Uik)9(p — k — pF)dk 
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Bi / G+(p — k,Uik)dk 
Jo 







A,/ G + (p-k,~ Ul k)9(p F -p + k)dk} (10) 



In the limit p ^ pp the last two expressions are equal and the limit is 



1 /" 

lim D(p) = — Sj[ / G + (p F - k,Uik)dk 
P->Pf Ztt ' J 

t—c.s ^ 

poo 

— / G+(pf — k, Uik)dk] = — 1 
Jo 



(11) 



above and bellow pp. Comparison of the exact expression for the momentum distribution function and its expansion 
around pp (both shown later) shows that 

D(p = pp) = — 1 + const. \p — pp \ a +const.(p — pp) (12) 

Therefore 

G+(p^pp,0) ~ -const } P ~ PF ^ (13) 

and therefore the Green's function changes sign at p — pp which completes the proof of the first part of Luttinger's 
theorem. 

Next, we would like to show that the Green's function given by Eq.(Q) gives the same number of particles as 
the noninteracting one. We use a band width cutoff A. Practically, that means that when the Green's function is 
integrated over the momentum, one integrates in the interval [pp — A,pp + A] and then takes the limit A/pp — *■ 
where A = pp — A. The number of particles will be the same if and only if: 

1 n ,___ f f dpde „ _s U t 



***AJ Jj2^ Go+Mel£ x 

dkduG + (p-k,e-Lj)K{k,uj) = (14) 



since 



— (l-— ) = » + (15) 



n \ pp 

here n + is the density of particles moving to the right. 

Using the mentioned properties of the Green's function, first we integrate over the frequencies e. Because of the 
second property, the integrals on the l.h.s. of Eq.(^) are convergent and therefore one can introduce the limit under 
the integral and perform it explicitly before integrating. This allow us to close the contour of integration on either 
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side of the real frequency axis in the complex plane. Therefore, when the branch-cut singularity of Green's function 
and the pole occur on the same half plane we close the contour in the other half plane and from Cauchy theorem 
the corresponding integral is zero. When the pole and the branch-cut singularity of the Green's function are on the 
opposite side of the real axis we choose to close the contour of integration in the half complex plane where the pole is 
located and therefore we obtain the branch-cut part evaluated at that pole. Then we evaluate the second frequency 
integration in ( |l4| ) using the same procedure. After performing the double frequency integration, the Ai terms give 
zero both in the charge and spin sectors. This is the reason why the simple square root Green's function |l2| leads 
to the step-function momentum distribution function, the same as for the non-interacting case. However, the double 
frequency integration of the Bi terms is nonzero and the result is: 

n + (p > p F ) = Bj J G+(p - k,v F (p ~ Pf) + u t k)0{p F -p + k)9(k)dk 

i—c,s 

n + (p <p F ) = 1 + — B * J ' G+(p-k,v F ( P -p F ) + u i k)0(p-k-p F )9(-k)dk (16) 

i— c,s 

The exact value of the last two integrals cannot be obtained without an explicit expression for the function. However, 



from the exact solution j27| and from perturbative calculations ]12|,|21 28 1 , it is known that the momentum distribution 



function is continuous and without a jump at p = p F and for p ~ p F is 

n(p>pr)~--&\p- p F \ a -C 2 (p - p F ) 

n(p< PF )~ l - + C 1 \p- p F | Q -C 2 (p - PF ) (17) 

where C\ and C 2 are constants |^8| . From the last two representations of the momentum distribution function follows 
the behavior of the integrals of the type encountered in the expression for G+(p, 0) in the neighborhood of p F . 

In the perturbation theory accessible regime (a < 1) the derivative ^ ~ \p — p F \ a ~ x is singular approaching infinity 
with a power law. We adopt the general definition that the Fermi surface is the set of /c-points at which the mth 
derivative of the momentum distribution function has a singularity i.e., {p F } = {Vfc : d d ^m^ is singular}. These we 
shall call Fermi points of order m. In the usual Fermi liquid the Fermi surface consist of zero order Fermi points while 
in the Luttinger liquid these are of first order. Although the zero order Fermi surface has disappeared, the first order 
Fermi surface is left and the generalized statement of the second part of Luttinger's theorem holds. When a > 1 
liquid droplets form and as long as a is not an integer the derivative is zero at p Fj but there will exist a number m 
so that is singular at p F and this will correspond to a Fermi surface of order m. 

At the end we must show that the following expression is zero in the limit A — > 0: 

p F p-p f 

I= ~^Y2 B ^J dp J d qG+(q,v F (p- p F ) +u l (p~ q)) 



t—c,s 



A 



+ J dp J dqG + (q,v F (p-p F )+Ui(p-q))\ (18) 
Pf p-p f 

Taking into account the two representations of the momentum distribution function one sees that in the above 
mentioned limit the two integrals cancel each other and the total number of particles on the right branch is n + = p f /tt. 

In this paper, we have shown from Dyson's equation, that Luttinger's theorem holds for the one-dimensional 
Tomonaga-Luttinger model. In general, the theorem is based on the counting of the fermionic degrees of freedom 
before and after the interactions are turned on. In the case of Fermi liquids, the one-to-one correspondence between 
the noninteracting particles and the quasiparticles ensures the validity of the theorem. In the case of ID Luttinger 
liquids, the number of charge particles in the interacting system is exactly equal to the number of electrons in the 
noninteracting system. Our conjecture is that as long as the number of states and excitations with and without 
interactions are the same, Luttinger's theorem will be satisfied. Some recently proposed 2D Luttinger-type Green's 
functions satisfy this condition The HTSC's, at zero doping and bellow the Neel temperature are antiferromag- 

netic insulators. Doping destroys the antiferromagnetic order and a metallic phase occurs above the superconducting 
critical temperature. Our discussion indicates that in the high-temperature superconducting cuprates, Luttinger's 
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theorem will be satisfied in this temperature and doping interval regardless of the nature of the electronic liquid 
(Fermi or Luttinger). 

We would like to thank J. Voit for the detailed analysis of this paper and the helpful suggestions. Special thanks 
are also due to J. Engclbrecht and M. Gulacsi for the insightful suggestions and to A. Nazarenko for carefully reading 
the paper. We would like to thank P.W. Anderson, C. DeLeonc, I.E. Dzyaloshinskii, and G. Zimanyi, for the useful 
discussions of this topic. 
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